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ABSTRACT 

The first low radio frequency (< 1.4 GHz) detection of the outburst of the 
recurrent nova RS Ophiuchi is presented in this letter. Radio emission was de- 
tected at 0.61 GHz on day 20 with a flux density of ~ 48 mJy and at 0.325 GHz 
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on day 38 with a flux density of ~ 44 mJy This is in contrast with the 1985 
outburst when it was not detected at 0.327 GHz even on day 66. The emission 
at low radio frequencies is clearly non-thermal and is well-explained by a syn- 
chrotron spectrum of index a ~ —0.8 (S oc u a ) suffering foreground absorption 
due to the pre-existing, ionized, warm, clumpy red giant wind. The absence of 
low frequency radio emission in 1985 and the earlier turn-on of the radio flux in 
the current outburst are interpreted as being due to higher foreground absorption 
in 1985 compared to that in 2006, suggesting that the overlying wind densities 
in 2006 are only ~ 30% of those in 1985. 

Subject headings: binaries: close - novae, cataclysmic variables - stars: individual 
(RS Ophiuchi) - stars: supernovae - stars: radio continuum 



1. Introduction 



The recurrent nova R S Ophiuchi was discovered to be in outburst on 2006 February 12.83 
UT (jNarumi et al. 20061 ). reaching a magnitude o f V = 4.5. The previous recorded out bursts 
took place in 1898, 1933, 1958, 1967 and 1985 faosino 19871: iRosino & Iiiima 1987h . with 
possible outbursts in 1907 JSchaeffer 2004h and 1945 ffoppenheimer & Mattei 1993h . 



The interacting binary system of RS Oph comprises an 



VI giant and a hot accret- 



ing white dwarf with an orbital period of 455.72 ± 0.83 days (IDobrzycka fc Kenyon 1994 



Anupama &: Mikolajewska 19991 : iFekel et al. 20001 ). The nova outbursts are powered by a 
thermonucle ar runaway on the white dwarf surface following accretion of mass from the 
companion fjStarrfield et al. 19851 : iKato 19901 ). There is remarkable similarity between the 
optical light curves and spectra from different outbursts (e.g. Rosino 1987) . The 1985 
outburst of RS Oph was one of t he best stud ied events, with the outburst being recorded 
from X-rays to radio wavelengths (IBode 19871 ). 



The first detection of the radio emission from RS Oph was made 18 days after the 1985 
outburst on January 26, by Padin et al. (1985) who observed a rise in the flux density of the 
radio source over ~ 20 days and noted that the data indicated high brightness temperatures 
which suggested a nonthermal origin. Monitoring of the 1985 outburst wi th the VLA at 
frequencies of 1.49, 4.85, 4.885, 14.94 and 22.46 GHz jHiellming et al. 1986h indicated that 
at least two radio components existed from roughly one to six months; one with a negative 
spectral index (where S u oc u a ) between 1.4 and 5 GHz; and another with a positive spectral 
index above 5 GHz. The positive index component completely dominated the decaying 
spectrum. Taylor et al. (1989) suggested the emission to be thermal at the higher frequencies. 
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An interesting feature to be noted in the RS Oph radio source is that while it was easily 
detected at all the high er frequencies (> 1.4 GHz) observed, it was not detected at 0.325 GHz 
( Spoelstra et al. 19871). The emission from the 1985 outburst was also extensively modelled 
f lBode & Kahn 19851 : lo' Brien et al. 1992h . 



The 2006 outburst of RS Oph has been subject to intense monitoring in all wavebands, 



Sokoloski et al. 2006 



almost immedia tely after discovery (10' Brien et al. 20061 ; iBode et al 2006 ; 
Das et al. 20061 ). Th e radio emission w as detected, at frequencies > 1.4 GHz as early as 4.7 
days from outburst ( Eyres et al. 20061). The radio source of RS Oph was resolved on day 
13.8 in the VLBI observations (10' Brien et al. 20061 ). The source, which was found to be an 
almost complete ring in the initial observations, soon evolved to a complex, multi-component 
structure consisting of an equatorial ring and polar caps. The early asymmetries in the radio 
source were explained by O'Brien et al. (2006) as being due to foreground absorption. Based 
on the observed brightness temperature and the estimate of the density in the shell, O'Brien 
et al. (2006) concluded that the radio emission is dominated by a non-thermal synchrotron 
component. 

We present in this letter the low frequency radio observations of RS Oph using the Giant 
Metrewave Radio Telescope (GMRT). We report the detection, for the first time, of RS Oph 
at frequencies below 1.4 GHz, and interpret the observed light curves. 



2. Observations 



Continuum observations with the GMRT (jSwarup et al. 19911 ) were first made on 2006 
February 24 (day 12.62) at 1.28 GHz. Subsequently, monitoring of the nova was performed at 
0.61, 0.325, 0.24, and 0.15 GHz. Since the 0.61 and 0.24 GHz feeds at GMRT are concentric, 
we observed both the frequency bands simultaneously. The last observation at 0.61 GHz was 
on 30 Jan 2007, and on 24 Nov 2006 at 0.325 GHz. 3C286 was used as the flux calibrator. 

The data obtained in native 'LTA' format were converted to standard FITS format and 
analysed using standard tasks in NRAO AIPS3. Self calibration using wide-field imaging was 
used to improve the image quality and, on the average, two rounds of phase self calibration 
were found to give the best images. Radio emission was detected at all the observed radio 
frequencies, except at 0.15 GHz observed on 2 June 2006 (day 110) to a 3a limit of 27 mJy 
We detected radio emission from the source at 0.325 and 0.61 GHz on all the observed days. 



1 AIPS is distributed by NRAO which is a facility of the NSF operated under cooperative agreement by 
Associated Universities, Inc. 
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The typical beamsizes at these two frequencies were ~ 10" and ~ 6" respectively and the 
source was not resolved in our observations. The uncertainty in the flux scale is ~ 15% and 
the errors on the flux densities listed in Table [T] reflect this uncertainty. 



3. The Radio Light Curves 

Fig. [1] shows the flux density evolution of RS Oph at 0.24, 0.325 and 0.61 GHz observed 
with the GMRT in the period 20 — 351 days after the outburst. The first observations at 0.61 
GHz, 0.325 GHz and 0.24 GHz were on March 5 (day 20), March 23 (day 38) and March 30 
(day 47) respectively. Also plotted in Fig. [T] are the L band observations using GMRT and 
the 1.49 GHz observations from the 1985 outburst. 

The light curves indicate a steep rise in the flux density, followed by a relatively flat 
maximum and a subsequent decay. While the different frequencies seem to become visible 
and peak at different epochs with the lower frequencies turning on at later times which is clear 
from our upper limit at 0.24 GHz around day 25 and subsequent detection around day 45; 
the post-maximum decay at all the frequencies is fairly similar. The observed spectral index 
varies from a ~ —0.1 around maximum to a ~ —1.0 around day 220. These observations 
clearly show that the radio emission at these frequencies is non-thermal. Although the non- 
thermal nature of the radio emission was inferred from the brightness temperatures and 



comparison between radio and X-ray flux (ITaylor et al. 19891 ) for the 1985 outburst, this 



is the first time that the non-thermal nature of the low frequency radio emission has been 
clearly demonstrated by the negative spectral index. 



4. Discussion 



The turn-on delay at longer wavelengths and the power-law decline after maximum 
with index (3 and a decreasing spectral index a observed in RS Ophiuchi are very simi- 
lar to the properties of radio supernovae, and imply the emission to be nonthermal syn- 
chrotron. We model the observed lig ht curves adopting the models for radio supernovae 
(jWeiler et al. 19861 : IWeiler et al. 20021 ) wherein the relativistic electrons and enhanced mag- 
netic fields necessary for synchrotron emission are generated due to the shock interaction of 
the nova ejecta with t he circumbinary red giant wind m aterial that is ionized and heated 
by the nova explosion (IChevalier 1982al ; IChevalier 1982bl ) . The radio emission rises rapidly 
as the shock progressively overtakes the wind material, causing a decrease in the line of 
sight absorption. Evidence for the presence of such an ionised, warm red giant wind is 
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provided by the X-r ay (IMason et al. 19871 ; iBode et al 20061 ; ISokoloski et al. 20061 ). radio 
( 10' Brien et al. 20061 ) and optical (G.C. Anupama et al. 2007, in preparation, Bode et al. 
2007) observations. 



4.1. Model light curves 

The following model, based on Weiler et al. (2002) is adopted. 



S(mJy) = K x 



/ _ r CSM 

V \ a I t — to\ CSM I 1 — e clumps 

g homog 



1GHz/ V 20d 
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T. 



CSM 
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(1) 

(2) 
(3) 



In the above model, to corresponds to 2006 Feb 12.83 and t — t is time since the 
outburst. Ki represents the flux density, K 2 the attenuation by a homogeneous absorbing 
medium and K 3 the attenuation by a clumpy/filamentary medium, at a frequency of 1 GHz, 
20 days after the nova explosion. The optical depths, T^^ g and t^^L are due to the ionized 
circumstellar material (CSM) external to the emitting region, a gives the spectral index and 
(3 the rate of decline in the optically thin phase. The optical depths in the homogeneous and 
clumpy /filamentary CSM are described by 5 and 5' respec tively. A simultan eous, non-linear 
chi-square fit using the Levenberg-Marquardt algorithm (IPress et al. 2002l ) is made to the 
data at all the frequencies and the best fit parameters (x re d = 1.5) are determined as listed 
in Tabled The flux density at 1.46 GHz on day 4.7 (2.8±0.2 mJy), reported by Eyres et al. 
(2006), is also used in the model fit. The resulting model light curves are plotted in Fig. [TJ 
It is clear that the overall evolution of the radio emission at freq < 1.4 GHz is well described 
by the model. The model predicts the delayed turn-on at different frequencies fairly well 
(ref . Table [2]) , implying the turn-on time is determined by the thermal optical depth of the 
CSM. However it underpredicts the peak fluxes at the 0.325 GHz and 0.24 GHz. Using our 
model fits and Eqns (11) and (13) of Weiler et al. (2002), w e estimate a mass lo ss rate of 
1.6 x 1O~ 7 M yr _1 for a red giant wind velocity of 20 kms -1 (IBode fc Kahn 19851 ). 



Excess emission, over the model fit, is observed at 0.61 and 0.325 GHz beyond day 
120. No such flux density variation is observed in two other bright sources in the field, 



- 6- 



implying the enhancement is intrinsic to RS Oph. It is interesting to note that S. P. S. Eyres 
et al. (2007, in preparation) record a slight enhancement in the flux density at 1.4 and 5 GHz 
on day 155. The enhancement that we record appears to be different from the short duration 
(< 24 hrs) increase in the flux density that was observed at 5 GHz, 41 days after the 198 5 
outburst and interpreted as being due to increased thermal emission (ISpoelstra et al. 19871 ). 
To quantify the difference in the observed and modelled light curves, we followed Ryder et 
al. (2004) and estimated the deviations for all the frequencies. Although the residuals are 
noisy, excess emission by a factor of 1.8 to 1.4 is observed between days 120 and 200. The 
deviations in the light curve may be interpreted as being due to modulations in the CSM an d 
the deviations are expected to be pronounced for edge-on viewing of the binary (IRvder 2004T) . 
The radio luminosity is related to the change in the CSM density (jChevalier 1982af ) and for 
the RS Oph system (with an inclination of 30° — 40°; Dobrzycka &; Kenyon, 1994 ), the 
observed enhancement can be explained if the CSM density is increased by 50% to 25%. 
It is interesting to note that although it has been proposed that the forward shock may 
traverse the red giant wind in around 80 days (e .g Bode et al. 2006 ), VLBI and HST 
observations (10' Brien et al. 20061 ; iBode et al 20071 ) are consistent with enhanced density in 
the pre-outburst CSM in the plane of the central binary. Shock breakout would therefore 
be expected to occur much later from material in this plane and thus this could be where 
the density enhancements inferred from our results reside. Alternately, the flux density 
enhancement could be related to the presence of a new short lived (~ 100 days) component 
in the radio emitting medium. The magnetic field for such a component, estimated assuming 
equipartition and the break frequency to be 0.325 GHz would be ~ 3.5 G. This is much 
higher than the equipartition magnetic fields in t he remnant which were estimated in th e 
1985 outburst to be of the order of 0.01 - 0.05 G faode & Kahn 19851 : iTavlor et al. 1989h . 
Spoelstra et al. (1987) have noted that a uniform magnetic field can be enhanced by a factor 
of 200 if the cooled remnant material is compressed in volume by a factor of 100 as the 
nova remnant expands. When the shocked medium becomes well cooled, Rayleigh- Taylor 
instabilities can cause such clumps to form. 



4.2. Comparison with 1985 radio light curves 



In Fig. [TJ the light curve for 1.49 GHz from the 1985 outburst (jHjellming et al. 19861 ) 
is plotted alongwith with the present low frequency light curves. The sharp rise to peak at 
0.325 and 0.24 GHz is very similar to the sharp rise to peak observed in 1985. The peak 
at 0.61 GHz is broader than at other frequencies. Although differences are seen around the 
maximum phase, it is interesting to note the very similar evolution during both outbursts in 
the post-maximum phase, particularly after day ~ 60. 
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We note two major differences in the radio emission during the 1985 and 2006 outburst: 



(i) The detection of emission at frequencies < 1.4 GHz. RS Oph was not detected at 
0.325 GHz on day 48 down to a la limit of 5 mJy during the 1985 outburst (jSpoelstra et al. 19871 ). 
while it was detected close to maximum on day 20 at 0.61 GHz and on day 38 at 0.325 GHz 
during the current outburst. 



(ii) The detection of radio emission as early as day 4.7 (lEyres et al. 20061 ) compared to 
an implied turn-on at day 14 in 1985 (e.g. Padin et al. 1985). 

Both these observations suggest that a steep-spectrum synchrotron source became visi- 
ble at the lower frequencies during the current outburst due to reduced foreground absorption 
compared to the 1985 outburst. Using the 3cr limit of 12 mJy on day 56 for the 0.325 GHz 
emission in 1985 (jSpoelstra et al. 19871 ). the observed value of 57 mJy on day 53 in 2006 and 
assuming that only foreground uniform absorption is responsible for the difference, we find a 
relation for the optical depths at 0.325 GHz in the two epochs r(1985) > r(2006) + 1.6 from 
simple radiative transfer arguments. The light curve model indicates that that r(2006) ~ 0.2 
on day 53 which implies r(1985) ~ 1.8. Since r oc n 2 e L, this indicates that the electron 
density of the foreground gas in 2006 is about 30% of that in 1985 for the sam e linear 
depth L. Since the two outbursts are otherwise fairly similar (jO' Brien et al. 20061 ). one of 
the reasons for this difference could be the variation in the dumpiness of the foreground 
gas. Observations of RS Oph at quiescence do indicate the presence of non-uniform col- 
umn density in the red giant win d envelope, that is uncorrelated with the binary geometry 
( lAnupama fc Mikolajewska 19991 ) . The non-detection at 0.325 GHz of the 1985 outburst had 
been attributed to absorption by the thermal gas mi xed with the synchrot ron emitting gas 
and low frequency cut-off in the electron spectrum (jSpoelstra et al. 19871 ). However both 
the physical processes appear to be negligible in the 2006 outburst. 



5. Summary 

GMRT observations of RS Oph during the 2006 outburst have shown early emission at 
frequencies below 1.4 GHz for the first time. The emission is clearly non-thermal, synchrotron 
emission and the evolution appears to be similar to that observed for radio supernovae. The 
light curves at these low frequencies are well explained by decreasing free-free absorption by 
the foreground CSM of the synchrotron emission from the remnant of the nova outburst. 
Model light curves with a spectral index of a = —0.8, decay power law of index (3 = —1.2 
and optical depth due to a homogeneous, clumpy absorbing medium well fit the observed 
light curves. Clumpy medium seems to dominate the absorption at early times. The light 
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curves indicate th e appearance of emiss ion components consistent with the VLBA, EVN and 



MERLIN images (10' Brien et al. 20061 ) at higher frequencies. Modelling this data with data 



at higher frequencies will enhance our understanding of this unique recurrent nova system. 

The similarity of late time evolution of the nova during the 1985 and 2006 outbursts 
suggests that the early difference is primarily due to the free-free absorption in the ionized 
CSM. This makes a case for early multifrequency observations with good temporal sampling 
during future outbursts. 



We thank the anonymous referee for insightful comments on the light curve modelling 
which have enhanced and improved the results of this paper. We thank Prof. R. Nityananda 
for providing TOO time on the GMRT during the early stages of the nova outburst and NGK 
thanks him for discussions. We thank the staff of the GMRT that made these observations 
possible. GMRT is run by the NCRA, a centre of the TIFR. This research has made use of 
NASA's Astrophysics Data System. 
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Table 1. Observed flux densities for RS Ophiuchi. The first column indicates the day 
after the outburst, 2006 Feb 12.83 UT (JD 2453779.33). 



Day 


V 


S 




Phase cal 




GHz 


mJy 


Name 


S (Jy) 


17.19 


1.39 


56.8 ±8.5 


1743- 


-038 


3.5 ± 0.5 


11.22 


1.28 


49.5 ±7.4 


1743- 


-038 


2.7 ± 0.4 


17.19 


1.28 


50.0 ± 7.5 


1743- 


-038 


O T _l_ f\ A 

2.7 ± 0.4 


17.19 


1.06 


55.4 ±8.3 


1743- 


-038 


2.6 ± U.4 


20.26 


0.61 


48.4 ± 7.2 


1743- 


-038 


1.1 ± U.Z 


29.15 


0.61 


48.9 ± 7.3 


1833- 


-210 


9.8 ± 1.5 


45.07 


0.61 


47.9 ± 7.2 


1822- 


-096 


7.7 ± 1.2 


92.99 


0.61 


21.0 ±3.2 


1822- 


-096 


8.1 ± 1.2 


i on sfi 


U.Ol 


10. It Z.o 


1822- 


-096 


O A _L 1 O 

8.4 ± 1.3 


136.92 


0.61 


15.3 ±2.3 


1822- 


-096 


or i t o 

8.5 ± 1.3 


147.88 


0.61 


10.5 ± 1.6 


1822- 


-096 


O O 1 T O 

8.8 ± 1.3 


153.94 


0.61 


18.2 ±2.7 


1822- 


-096 


9.0 ± 1.4 


166.77 


0.61 


10.9 ± 1.6 


1822- 


-096 


inn _i_ i c 
10.0 ± 1.5 


192.84 


0.61 


10.1 ± 1.5 


1822- 


-096 


9.9 ± 1.5 


206.89 


0.61 


5.5 ±0.8 


1822- 


-096 


8.2 ± 1.2 


214.66 


0.61 


5.8 ±0.9 


1822- 


-096 


8.3 ± 1.2 


ZZl.oO 


U.Ol 


4.Z Zt Z.l 


1822- 


-096 


9.2 ± 1.4 


351.40 


0.61 


3.9 ±0.6 


1822- 


-096 


7.8 ± 1.2 


38.15 


0.325 


43.7 ±6.6 


1822- 


-096 


10.4 ±1.6 


53.24 


0.325 


57.0 ±8.6 


1822- 


-096 


11.8 ±1.8 


123.12 


0.325 


17.2 ±3.4 


1822- 


-096 


9.8 ±1.5 


134.07 


0.325 


18.0 ±3.6 


1822- 


-096 


10.7 ± 1.6 


169.75 


0.325 


20.6 ±5.7 


1822- 


-096 


10.1 ± 1.5 


225.61 


0.325 


8.1 ±4.4 


1822- 


-096 


9.5 ± 1.4 


284.43 


0.325 


6.3 ±1.8 


1822- 


-096 


9.7 ± 1.5 


29.15 


0.24 


< 13.0 








45.07 


0.24 


54.2 ±8.0 


1822- 


-096 


14.0 ±2.1 
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Table 1 — Continued 



Day 


V 


S 


Phase cal 




GHz 


mJy 


Name S (Jy) 


92.99 


0.24 


30.5 ±4.6 


1822-096 15.0 ±2.3 


110 


0.15 


< 9.0 





Table 2. Best fit parameters from least square fits to the low frequency radio emission 

from RS Ophiuchi 



Parameter 


Value 




86.4 


a 


-0.78 


P 


-1.24 


K 2 


0.144 


5 


-2.29 


K 3 


0.53 


5' 


-3.14 


v (GHz) 


Onset (days) 


1.28 


3 


0.61 


6 


0.325 


10 


0.24 


14 



-14- 




log(t-t ) day 



Fig. 1.— Observed light curves of RS Oph (points) at L-band, 0.61 GHz, 0.325 and 0.24 GHz 
observed in the 2006 outburst. Continuous lines represent the model light curves generated 
using par ameters listed in Tab le [21 Also plotted is the 1.49 GHz light curve from 1985 
outburst (iHjellming et al. 19861 ). 



